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Abstract 

We study critical spreading dynamics in the two-dimensional contact process 
(CP) with quenched disorder in the form of random dilution. In the pure 
model, spreading from a single particle at the critical point A c is characterized 
by the critical exponents of directed percolation: in 2+1 dimensions, 5 = 0.46, 
T) = 0.214, and z = 1.13. Disorder causes a dramatic change in the critical 
exponents, to 5 ~ 0.60, r/ ~ —0.42, and z ~ 0.24. These exponents govern 
spreading following a long crossover period. The usual hyperscaling relation, 
45 + 2r] = dz, is violated. Our results support the conjecture by Bramson, 
Durrett, and Schonmann [Ann. Prob. 19, 960 (1991)], that in two or more 
dimensions the disordered CP has only a single phase transition. 

PACS numbers: 05.50.+q, 02.50.-r, 05.70.Ln 
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Phase transitions between an absorbing state (one, that is, admitting no further evolu- 
tion), and an active regime occur in models of autocatalytic chemical reactions, epidemics, 
and transport in disordered media. Critical phenomena attending absorbing-state transi- 
tions show a high degree of universality, that has been characterized rather precisely in 
studies of the contact process (CP) and of directed percolation (DP) [T|-§[. Since many- 
particle systems often incorporate frozen-in randomness, it is natural to investigate its effect 
on an absorbing-state transition. Some time ago, Noest studied the change in the static 
critical behavior of DP due to quenched disorder Q. In this work we examine its effect 
on time- dependent critical phenomena at an absorbing-state transition. We focus on the 
two-dimensional CP, a simple lattice model of an epidemic ||. Our prime interest is in the 
effect of disorder on the spread of the critical process from a seed. 

At a critical creation rate, A c , the pure CP exhibits a second-order phase transition 
characterized by the same critical exponents as DP [0J . The well-known Harris criterion ||[7] 
states that disorder changes the critical exponents if du± < 2, where d is the dimensionality 
and u± is the correlation-length exponent of the pure model. Since v± ~ 0.73 for DP in 
2+1 dimensions, we expect quenched disorder to be relevant in the CP. Indeed, Noest 's 
simulations of one- and two-dimensional stochastic cellular automata (belonging to the DP 
class), yielded static critical exponents quite distinct from those of DP, when the models 
were modified to incorporate quenched randomness ||. A field-theoretic study by Obukhov 
H yielded qualitatively consistent results. Marques studied the effects of dilution on the CP 
and related models in a mean-field renormalization group study, obtaining a phase diagram 
in good agreement with simulation |J. 

The effect of disorder on time- dependent critical behavior at an absorbing-state transition 
has received little attention. The only study we are aware of, by Bramson, Durrett, and 
Schonmann, showed that the one-dimensional CP in a random environment has a new phase 
in which the process survives but spreads more slowly than linearly with time, t |I0| . (In 
the pure CP the radius of the active region grows oc t for any A > A c .) 

In this paper we report on extensive simulations of time-dependent critical behavior in 
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the two-dimensional DCP. The model is defined as follows. Each site of the square lattice 
Z 2 is either vacant or occupied by a particle. Particles are created at vacant sites at rate 
An/4, where n is the number of occupied nearest-neighbors, and are annihilated at unit 
rate, independent of the surrounding configuration. The order parameter is the particle 
density p; it vanishes in the vacuum state, which is absorbing. As A is increased beyond 
A c = 1.6488(1), there is a continuous phase transition from the vacuum to an active steady 
state. We introduce disorder by randomly removing a fraction x of the sites. That is, for 
each G Z 2 there is an independent random variable taking values and 1 with 

probability x and 1 — x, respectively. The DCP is simply the contact process restricted 
to sites with r)(i,j) = 1; those having f](i,j) = are never occupied. (Thus if exactly m 
neighbors of a given site have r)(i,j) = 1, the creation rate at that site is at most mA/4.) 
Naturally, 1 — x must exceed the percolation threshold p c = 0.5927 for there to be any 
possibility of an active state, since on finite sets the CP is doomed to extinction. 

Following the approach of Grassberger and de la Torre [I], we study a large ensemble of 
trials, all starting from a configuration very close to the absorbing state: a single particle at 
the origin. For A > A c (x) there is a nonzero probability that the process survives as t — > oo; 
for A < A c (x) the process dies with probability 1. Of prime interest are P{t), the survival 
probability at time t, n(t), the mean number of particles (averaged over all trials, including 
those that die before time t), and R 2 (t), the mean-square distance of particles from the 
origin. At the critical point of the pure CP, these quantities follow asymptotic power-laws, 

P(t)oct- 5 (1) 
n(t) cx f (2) 
R 2 (t) cx t z . (3) 

The exponents satisfy the hyperscaling relation 

A5 + 2rj = dz, (4) 

in d < 4 dimensions [||]. For A < A c , P(t) and n(t) decay exponentially, while for A > A c , 
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P(t) -»• Poo ~ (A - X c f , (5) 

and n(t) ~ t d . Thus log-log plots of P(t),n(t) and R 2 {t) approach straight lines at the 
critical point, and show a positive or negative curvature in the supercritical or subcritical 
regimes, respectively. We adopt the working hypothesis that the critical point of the DCP 
is also characterized by asymptotic power laws. 

We studied dilutions x = 0.02, 0.05, 0.1, 0.2, 0.3, and 0.35. Samples of from 10 5 to 2x 10 6 
trials were generated for each A value of interest, each trial extending to a maximum time 
tmax < 7 x 10 4 . (As is usual in this sort of simulation, the time increment associated with 
an elementary event — creation or annihilation — is At = 1/N, where N is the number 
of particles.) An independent realization of disorder (the variables 7](i,j)), is generated 
for each trial. Our results for P(t), n(t) and R 2 (t) (see Fig. 1), allow us to compute the 
local slopes 8(t), r](t), and z(t), using least-squares fits to the data (in a logarithmic plot), 
distributed symmetrically about a given t (typically in the interval [t/3, 3t]). We estimate 



as 



the exponents by plotting the local slopes versus 1/t and extrapolating to 1/t — ► [11 
in Fig. 4. The critical point A c (x) is determined by the criteria that at long times, log-log 
plots of P(t) appear straight, and that local-slope plots are free of marked curvature. Our 
results for the critical point are given in Table I. For small x, X c {x) ~ A c (0)/(1 — x), as 



predicted by mean-field theory P,[12 



Perhaps the most surprising feature of the critical DCP is the population size n(t). In 
the (pure) critical CP n(t) increases monotonically with time; rj in Eq. (2) is positive. In the 
critical DCP n(t) decays at long times; rj is negative. The critical point evolution of P, n, 
and R 2 , shown for several dilutions in Fig. 1, has a regime characterized by pure CP scaling, 
followed by a long crossover period leading to a different set of power laws. The apparently 
disparate results obtained for different dilutions are in fact described by a single function. 
For n(t), this can be seen by shifting the data (plotted on log-log scales), to superimpose 
the maxima, which causes all of the points to fall close to a single curve (see Fig. 2). If r(x) 
denotes the time at which n(i; x) attains its maximum, n max (x), then we have 
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n(t; x) = n max (x)M(t) , (6) 

where t — tjr and M is a scaling function. We regard r as the crossover time separating pure 
CP scaling (for small t) from disordered CP critical behavior (t — > oo). To scale the data for 
P(t) and R 2 (t) we use the same t(x) values as for n, and fix dilution-dependent amplitudes 
Ap(x) and A R 2(x) by optimizing the data collapse. Thus for the survival probability 

P(t;x) = A P (x)V(t), (7) 

and for the mean-square spread, 

R 2 (t;x) = A R2 (x)K{t) . (8) 

To test these scaling expressions we plot P = P(t;x)/A P (x), h = n(t; x)/n max (x) and 
R 2 = R 2 (t; x)/A R2 {x) versus t in Fig. 2. Data for different dilutions collapse nicely onto a 
single curve. In Fig. 3 we plot the crossover time r and the amplitudes versus x. For x > 0.1 
these parameters are well-described by power laws, for example, r cx x' 3 - 2 ^. (Figures in 
parentheses denote uncertainties. The corresponding exponents for Ap, n max , and A R 2 are 
1.67(4), -0.56(2) and -3.64(3), respectively.) Given the limited range of x, there is no reason 
to attribute any fundamental significance to these power laws. The rapid decrease in r 
with increasing x can be understood by noting that for weak disorder, the process must 
spread over a certain area before randomness becomes manifest; for small x, sizable regions 
of the lattice look nearly regular. Thus for x < 0.2, n(t) and R 2 (t) have not reached their 
asymptotic behavior at the end of the trials (t max = 7 x 10 4 ). Since n(t; x = 0.02) does not 
exhibit a maximum on simulation time scales, the data for x = 0.02 are not included in the 
scaling analysis. 

Given that the approach to asymptotic scaling is more rapid for larger dilutions, we turn 
first to the data for x = 0.35 to estimate critical exponents. Fig. 4, for r/(t), shows that 
the local slope undergoes a marked change during the evolution. (z(t) also changes over a 
very wide range.) At long times the local slopes for x = 0.3 and x = 0.35 exhibit steady 
trends, permitting reliable extrapolation of their asymptotic values; for smaller dilutions 
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only 5 can be estimated with confidence (see Table I). The scaling plots of Fig. 2 show 
that the scaling functions V, J\f, and 1Z follow power laws with pure-DP exponents for small 
t = t/r, and that for large t they approach power laws with the exponents found in the 
local-slope analyses. 

We also determined the ultimate survival probability exponent (3' (defined in Eq. (|5[)), for 
dilutions 0.1 and 0.2. While f3' is conceptually distinct from the order-parameter exponent 
P, defined through p ~ (A — A c ) , these exponents are equal in the CP and allied models. 
Stationary-state simulations support the relation (3 = (3' for the DCP as well ||12|| . 

Table I reveals the dramatic effect dilution has on the exponents. Comparing the critical 
exponents for the DCP with those of the pure model, it is evident that the shifts all reflect a 
diminished ability of the critical process to survive and expand in the presence of disorder. 
While this accords with intuition, the small z value, and negative r\ are unexpected. Most 
absorbing-state transitions have r\ > 0, but this does not appear to be a fast rule, so long 
as r] s = i] + 5 > 0. (r] s governs the population in surviving trials.) Negative values of 
7] were also found in a (disorder-free) model with multiple absorbing configurations ||T3| . 



Equally puzzling is the massive violation of hyperscaling, Eq. (|J). In fact, hyperscaling 
in the DCP may be incompatible with the simple requirement r] s > 0; for d — 2, Eq. (|J) 
implies r] s = z — 5, which in turn yields (for x = 0.35), r] s ~ —0.35. Alternatively, if we take 
1] = —5, (that is, as small as possible without forcing rj s < 0), Eq @ yields z = S, in marked 
disagreement with the data. 

One might regard violation of hyperscaling as evidence that our simulations have yet to 
probe the asymptotic scaling regime. While simulations which run to a finite time cannot 
rule this out, the data suggest no trend that would bring the exponents into agreement with 



hyperscaling. Another interpretation is that sample-to-sample fluctuations and lack of 
self-averaging lead to a violation of the scaling hypothesis itself (as set out, e.g., in Ref. JT|). 
We hope to resolve this issue in a future study. 

Bramson, Durrett and Schonmann studied a one-dimensional CP with disorder in the 



form of a death rate randomly taking one of two values (independently) at each site |T0 



They demonstrated that this model possesses an intermediate phase in which survival (start- 
ing, e.g., from a single particle) is possible, but the active region grows more slowly than 
linearly. In two or more dimensions, Bramson et al. conjectured, there is no intermediate 
phase. Our results support this conjecture. Simulations at x = 0.1, with A close to, but 
slightly above A c (to be precise, A = 1.86 and 1.87, corresponding to (A — A c )/A c = 0.008 and 
0.014, respectively), showed n(t) ~ t 2 (and similarly for R 2 (t)), consistent with the radius of 
the active region growing ~ t. Thus a sublinear growth phase, if it exists at all, is confined 
to a very narrow range of creation rates. While our model incorporates dilution rather than 
a random death rate, one would expect such an intermediate phase to be a rather general 
feature of disordered contact processes, so that its apparent absence here argues for the 
validity of the conjecture. 
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TABLES 

TABLE I. Critical parameters and exponents obtained through simulations of the DCP. The num- 
bers in brackets indicate uncertainties. 



X 


A c 


& 


5 


V 


z 





1.6488[1] 


0.586[14] 


0.460[6] 


0.214[8] 


1.134[4] 


0.02 


1.6850[3] 


0.566[7] 


0.467[1] 


0.216[3] 


1.104 [2] 


0.05 


1.7405 [3] 




0.59[1] 






0.1 


1.8448 [4] 


0.70[1] 


0.59[1] 






0.2 


2.097[2] 


0.72 [2] 


0.60[2] 






0.3 


2.437[2] 




0.59[1] 


-0.43[3] 


0.27[3] 


0.35 


2.655[2] 




0.59[1] 


-0.41 [3] 


0.24 [3] 
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In summary, we find that quenched disorder induces a radical change in the critical 
exponents describing spreading in the contact process. As the population spreads in the 
critical process, there is a crossover from pure CP behavior to a new set of critical exponents. 
The crossover time appears to diverge as the dilution tends to zero. A related finding is the 
demonstration by Noest of anomalously slow dynamics, in the form of power-law relaxation 
to the vacuum, in disordered directed percolation fl5fl . Although our results on scaling 
are restricted to dilutions 0.05 < x < 0.35, we expect that such a crossover occurs for all 
< x < 1 — p c , albeit at very long times for small x. One is inclined to suppose that the 
DCP is but one member of a universality class encompassing all disordered models with a 
continuous transition to a unique absorbing configuration. In support of this we observe that 
5 (the only exponent we are able to determine over a wide range of x values), is insensitive 
to the degree of dilution, and that evolutions for various x, with very different crossover 
times, are described by common scaling functions. Studies of absorbing state transitions in 
other disordered models are needed to verify the universality hypothesis. Finally, we note 
that our result, (3 — (3' — 0.71(2), is incompatible with that reported by Noest (/3 = 1.10(5)) 
for disordered DP in 2+1 dimensions [[|. We shall return to this question in a future work 
on the static behavior of the diluted contact process. 
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Figure Captions 

FIG. 1. Survival probability P, mean population n, and mean-square distance of particles 
from the origin R 2 , in the critical contact process with dilution x. o: x = 0.05; +: x = 0.2; 
x: x = 0.35. 

FIG. 2. Scaling plots of spreading in the DCP. Upper panel: InP vs lnt; the straight lines 
have slopes of -0.46 and -0.60. Middle: Inn; the slopes are 0.241 and -0.42. Bottom: lni? 2 ; 
the slopes are 1.134 and 0.24. Symbols: o: x = 0.05; □: x = 0.10; +: x = 0.20; A: x = 0.30; 
x: x = 0.35. 

FIG. 3. Log-log plots of r (o), n max (□), A P (+), and A R 2 (x) vs dilution x. 

FIG. 4. Local slope r)(t) versus t^ 1 in the diluted contact process with x = 0.35. Symbols: 
+ : A = 2.650; x: A = 2.655; o: A = 2.660. 
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